THIN DOMAINS WITH EXTREMELY HIGH OSCILLATORY BOUNDARIES 
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Abstract. In this paper we analyze the behavior of solutions of the Neumann problem posed in a thin 

domain of the type R"^ = {{xi,X2) G | xi G (0,1), — cb(xi) < x-z < cG{xi, xi/e")} with o > 1 and 
e > 0, defined by smooth functions b(x) and G{x, y), where the function G is supposed to be /(x)-periodic 
in the second variable y. The condition a > 1 implies that the upper boundary of this thin domain presents 
a very high oscillatory behavior. Indeed, we have that the order of its oscillations is larger than the order 
of the amplitude and height of R' given by the small parameter e. We also consider more general and 
complicated geometries for thin domains which are not given as the graph of certain smooth functions, but 
rather more comb-like domains. 



1. Introduction 



In this paper, we analyze the behavior of the solutions of the Laplace equation with homogeneous Neumann 
boundary conditions 

-Aw" = h" in R" 

where A'"'^ is the unit outward normal to dR'^ and h"^ € L'^{R"). The domain R" is a two dimensional thin 
domain which presents a highly oscillatory behavior at the boundary. We will be able to consider two 
different types of thin domains, which will be clearly defined in Section 2. To make the ideas clear we will 
refer in this introduction to the first type: assume ii^ is given as the region between two functions, that is, 

R" = {{xi,X2) e I a;i e (0,1), -eh{xi) <X2< eG,(a;i)} (1.2) 

where 6(-) and G'c(-) arc functions satisfying < 6o < ^(0 < < G^{ ) < Gi for some fixed positive 
constants bo, hi and Gi, independent of e > 0. Here, the function 6, independent of e, defines the lower 
boundary of the thin domain, and the function G^, dependent of e, the upper boundary of R'^. We will allow 
Ge to present oscillations whose amplitude is larger than the order of compression of the thin domain. This 
is expressed by assuming that 

G,(x) = G(a;,:c/e"), (1.3) 
for some positive constant a > \. The function G : (0, 1) x M — )• R is a positive smooth function, with 
y — > G{x,y) periodic in y for fixed x with period l{x). 

Let us observe that our assumptions includes the case where the function Gg presents a purely periodic 
behavior, for instance, G, (x) —2 + sin(x/e"). But it also considers the case where the function G^ defines 
a thin domain where the oscillations period, the amplitude and the profile vary with respect to a; € (0, 1). 
The Figure 1 and 2 below illustrate kinds of thin domains that we are considering here. 

Since the domain i?"^ is thin, approaching the interval (0, 1), it is reasonable to expect that the family of 
solutions will converge to a function of just one variable and that this function will satisfy certain elliptic 
equation in one dimension with some boundary conditions. 
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Figure 1. A thin domain with variable period, ampHtude and profile. 




Figure 2. A comb- like thin domain. 



It is known that if the domain does not present oscillations, that is Ge{x) — G{x), with < Go < G(-) < Gi 
and b{-) = the 1-dimensional limiting problem is given by 



:;^(G{x)Wa;) +W = h, ill (0,1), 
rix) V /x 



G(x)V"^-^'-^V. (1.4) 
u;,(0) = = 

see for instance (TUl [H]. Also, if we consider 6 = 0, G^{x) = G{x,x/e°') for some < a < 1, and if we 
assume that Ge(-) —> m(-) w — L^{0, 1) and k{-) w — L^{0, 1) (observe that m{x)k{x) > 1 a.e. 

and in general it is not true that m{x)k{x) = 1), then the limit problem is 

Wx\ + w ^ h, in (0, 1) 



m(x) \k{x) 

Wx{0) = wx{l) ^ 



see [2] for details. Note that this case contains the previous one since we can recover the problem (1.4) 
assuming a = 0. 

Recently, we consider in [H [5] a class of oscillating thin domain that cover the case a = 1 with constant 
period I. Observe that this situation is very resonant since the height of the domain, the amplitude of the 
oscillations at the boundary and the period of the oscillations are of the same order e. The limit problem 
for this case is 

-—-{r{x)wx)x+w = h{x), a;e(0, 1) 

s{x) (1.5) 
w'(0) = u;'(l) = 



where 

r{x) — 

IY'{x) 



|l- {yi,y2)\dyidy2 
Jy'(x) oyi J 

s{x) = \Y*{x)\ 
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and X{x) is a convenient auxiliary harmonic function defined in the representative basic cell y*(x), which 
depends on G(x, •), x £ (0, f ), and it is given by 

Y*{x) = {(yi,2/2) eR^\0<yi<l, < 2/2 < G{x,yi)}. 

The restricted case where the function Gi{x) — G{x/e) for some /-periodic smooth function G can be 
addressed by somehow standard techniques in homogenization theory, as developed in [51 [51 US] ■ We refer to 
[3] for a complete analysis of this case for a semilinear parabolic problem. 



In this work, we are interested in addressing the case a > 1 in (1.3 1, where none of the techniques used 
to solve the previous ones apply. In particular, we do not use any extension operator for the convergence 
proof. Indeed, we will be able to show how the geometry of the boundary oscillations affect the limiting 
equation, see Theorem |2.1[ Theorem |2.5[ See also Corollary |2.3| for a very interesting interpretation of the 
limiting equation and to see how the geometry of the unit cell affect the limit equation in the case of periodic 
oscillations. 

In Section[2|we give precise definitions of the two types of thin domains we are considering. One of them is 
the one described in this introduction. The other type is a "comb-like" thin domain, which can be visualize 
in Figure 2. We also state clearly the two main results we prove. Theorem |2 . 1 j and Theorem |2.5[ 

The short Section [3| states a technical result which will be used later in the proof. 

In Section [4| we analyze the type of thin domains which are given as a region between two graphs as in 



(1.2) 



In Section [5| we analyze the other type of thin domains, that we have denoted as a "comb- like" thin 
domain. 

Wc also would like to observe that although we will treat the Neumann boundary condition problem, we 
may also impose different conditions in the lateral boundaries of the thin domain R'^ , while preserving the 
Neumann type boundary condition in the upper and lower boundary. Indeed, we may consider condi tion s 



of the Dirichlet type, w*^ — 0, or even Robin, -I- (3w'^ = in the lateral boundaries of the problem (|l.l|). 
The limit problem will preserve this boundary condition. 

2. Basic facts, notation and main results 

We will consider two different types of thin domains. One of them will be given as the region between the 
graphs of two functions and the other will consists of an autoreplicating structure with appropriate scaling 
rates which resembles a comb structure. We present now the main definitions, basic facts and results on 
both cases. 

Type I. Thin domain as the region between two graphs. Let us consider a one parameter family of 
functions G^ : (0, 1) — )■ [0,oo), e € (0,eo) for some eq > 0, and a function b : (0, 1) (0,cx)). We will assume 
the following hypotheses on functions b and G^: 

(HI) There exist two positive constants 69, bi such that < bg < b{x) < bi for all x E (0, 1) and the 
function b is piecewise C^. . 

(112) The functions Ge{ ) are of the type G'e(x) = G'(x,x/e"), with a > 1, where the function 

G:[0,l]xM ^[0,+c^) 

{x,y) -^G{x,y) ^^■'> 

is continuous in x, uniform in the second variable y, (that is, for each 77 > 0, there exists S > 
such that \G{x,y) — G{x\y)\ < t] for all x, x' € [0,1], \x — x'\ < (5, and y € M). Moreover, we 
suppose G{x,y) > is periodic in y, with a period l{x) that may depend on the first variable, that 
is, G{x,y + l{x)) = G{x, y). We also assume that /(•) is a continuous function with < i' < l{x) < L 
for all X G [0,1]. 

We consider the highly oscillating thin domain i?*^, which is given as the region between the graphs of the 
two functions eb{-) and eGe(-), that is 

R' = {ixi,X2) e I xi G (0,1), ~£b{xi) <X2< eG,{xi)} 

and we investigate the behavior of the solutions of ( |1.1[ ) as e — > 0. 
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The Figure 3 below, gives us an exemple of function G in a bounded open set. 

Since a > 1, we have that the upper boundary of this thin domain presents a extremely high oscillatory 
behavior. More precisely, the order of the oscillations is large than the order of the amplitude and height of 
the thin domain i?*^ with respect to the small parameter e. Also, we get more general perturbations assuming 
that the period I depends on variable x E (0, 1). 




Figure 3. Graph of function G in (0, 1) x (0, 3L). 



To study the convergence of the solutions of (1.1 ), we consider the equivalent linear elliptic problem 

(2.2) 



^x^_'■ e2 dxi' 



ax\ ax2 
where S satisfies 

lirilL^(ao<C (2.3) 

for some C > independent of e, and now, = [v\^ f|) is the outward unit normal to 911^, and C is 
a highly oscillating domain given by 

r!^ = {(xi,X2) e I XI G (0, 1), -b[xi) <X2< G,{xi)}. (2.4) 



Note that the equivalence between ( 1.1 ) and ( 2.2 ) is easily obtained by changing the scale of the thin domain 
i?' in the y-direction through the simple transformation {x, y) — (x, ey), (see [21 HH] for more details). Thus, 
we have a domain which is not thin anymore but presents very wild oscillatory behavior at the top boundary, 
although the presence of a high diffusion coefficient in front of the derivative with respect the second variable 
decreases the effect of the high oscillations. 

We also mention the works [TJ [71 IS] that analyse elliptic problems in domains related to il' but the fact 
that in our case we allow very high diffusion in the y-direction with distinct oscillations makes our analysis 
and results different from these other papers. 

Now we are in contidion to state our main result whose proof will be presented in section |4.2[ 



rG,(x) 

,L J KZ -U yiL J C^Ubt^yJ iCCi \\J ||L-^(£1*^) — ^ lulu ULC- JLLlit^bLUIL J J 

satisfies that /"^ ^ /, w-L'^{0, 1). Let u'^ be the unique solution of (2.2) and Gq be the function given by 



Theorem 2.1. Assume that £ L'^{il'^) satisfies \\f^\\L^(n^-) < C and the function f"^ (x) — j'^'l^^^^ f{x,y) dy 



Go(a:) = minG(a;,2/) > 0. (2.5) 



Then, if Uq{xi) is the unique weak solution of the Neumann problem 



^1 

^(b{x) + Go{x)j u.^{x)ip^{x) +p{x)u{x)ip{x)^dx ^ J f{x)ipdx, V(/9 G i?\0, 1) (2.6) 
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where p{x) is the function defined as follows: 

p{x) - / h{x) + G{x, y))dy = h{x) + — / G(x, y) dy, for all x G (0, 1), 
l{x) Jq \ ) l{x) Jo 

we have 

Wu" - uoWl^si') ^ 0. 

Moreover, if we denote by fig — {{xi,X2) S | xi S (0, 1), ~b{xi) < X2 < Go{xi)} C fie then, 

Ue ^ uq w — i7^(r2o) 



Remark 2.2. The functions b and Gq defined in (HI) and (2.5) respectively, are associated to the part of 
the domain fl'^ that does not oscillate as the parameter e goes to zero, and have an important role in the 
limit problem (2.6). Indeed, if we assume that the period, the amplitude and the profile of the domain are 
constant with respect to x £ (0, 1), we get the nice result announced below. 

Corollary 2.3. // we have G{x, y) = G{y) an L-periodic function with mixYy^^ G{y) = and b{x) = b a 
constant function, then, the homogenized limit is given by the equation with constant coefficients: 

i~du" + u = f, (0,1) 

\ w'(0) =u'(l) = ^ ■ ^ 

where the diffusion coefficient is given by 

b Lb Area of the non oscillating part of the unit cell 

d — 

Remark 2.4. In Corollary \2.3\ the constant Lb represents the area of the unit cell which correspond to the 
part which is non oscillating and Lb + G{y) dy represents the total area of the unit cell. 

Type II. Comb-like thin domain. We consider now another interesting type of thin domain. Consider 

R" = Int (Rl U R^) 

where 

i?l = {{xi,X2) I < xi < 1, ~eb{xi) <X2< 0}, 
with b given as in the previous case (see hypothesis (HI)), and 

where 

i?^ + = {(nie" + e°'xi,ex2) \ (xi, Xa) € Qo} 
where Qo C (0,i) x (0,G) is a fixed Lipschitz domain satisfying the following: 

(HQ) if Fq = dQo n {x2 = 0} and if we consider ei(Qo) the first eigenvalue of the operator — A in Qo with 
homogeneous Dirichlet boundary condition in Fq and homogeneous Neumann boundary condition 
in dQo \ Fo, then ei((5o) > 0. 
Observe that if Qo is connected and Fq 7^ then (HQ) is satisfied. But there are cases where Qo is 
disconnected ant still (HQ) holds. 

As we have done in the previous case, let us define il^ = {(2^1, X2) \ [xi, 6X2) G i?'} so that, 



= Int (r2_ u n%) , 

n_ = {{xi,X2) I < xi < l,-6(a;i) < X2 < 0}, 

where 

= {{nLe" + X2) | (xi, X2) S Qo}. 
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We also consider the equivalent linear elliptic problem 





1 d^u" 








1 du" 




dx2 



+ U'' ^ f in 

v^ = Q on dVt" 



(2.8) 



Under this conditions, we may get the following result. 



Theorem 2.5. Let u'^ he the unique solution of {2.2). Assume that G i (51^) satisfies ||/'^||l2(si') ^ C 
and the function /'(a;) = Is'{x) /"^(^j y) "^2/ satisfies that /, w-L^(0, 1) where S'^{x) — {y \ (x, j/) G £7*^}, 

t/iaf is, t/ie section of the domain fl'^ at the point x S (0, 1). 

Then, if uq{xi) is the unique weak solution of the Neumann problem 



^b{x) Ux{x) (px{x) + q{x) u{x) ip{x)^ 
where q{x) is the function given by 



dx 



f{x)(pdx, V^ei7^(0,l) 



q{x) 



L 



+ b{xi) Va;ie(0,l), 



we have 
Moreover, 



- Uo\\l2(Q') ^ 0. 



3. An important estimate 

In this section we show several basic estimates on the solutions of certain elliptics pde's posed in rectangles 
of the type 

Q. = {{x, y)eR^\ - < X < e", < 2/ < 1} 
with a > 1. As a matter of fact, for Mo(-) G e",e"), we define the function u'^{x,y) as the unique 

solution of 

d^u'' 1 d'^u" 



dx' 



dy' 



in 



u(a;, 0) = uo(x), on Fj, 
du 



(3.1) 



dv 



0, on dQ, \ r. 



e" < a; < e"}. 



where v is the outward unit normal to dQf and 

r, = {(x,0)eM2 

We have the following. 

Lemma 3.1. With the notation from above, if we denote by uq the average of uq in T^, that is 

1 



Mo 



2e" 



uq{x) dx 



then there exists a constant C, independent of e and uq, such that 



\u%x,y) - uofdx < Cexpj--^ 



|tio|lL2(-£=,e°) 



1 /.£ 



\u{x,y) - uof dxdy < Ce" ^||uo|li2(_e=_,c) 



(3.2) 

(3.3) 
(3.4) 
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and 



du 


2 


1 


9u 


dx 






9y 



9uo 



9a; 



(3.5) 



Proof. The proof of this result is based in the known fact that the sohition of the problem above can be 
found explicitely and admits a Fourier decomposition of the form 

cosh(^^^fe^^) 



1 
2^ 



{x)dx + ^(uo, <y9;)(/?;(x)- 



k=l 



cosh(-^) 



where (^^(x) e "/2cos(^p5), n = 1,2, . . . , and (uo,0 = ("o, ¥'^)L2(_,= ,e°)- 



Remark 3.2. Observe that in particular, estimate (3.5) implies that 



min < 


9m 


1 


9u 


1 < Ce"-i 


9uo 


wev 1 






9?/ 




dx 


w/iere = {u G H^iQt) u = 


= Wo 


in Fe}. 









(3.6) 

□ 



4. Thin domains as a region between graphs. 

In this section we consider Type I thin domains and provide a proof of Theorem |2.1[ 
We will start analyzing in detail the structure of the domain fi*^ as a preparation for the proof of our 
result. 

4.1. The one parameter family G^. In this subsection we obtain some properties and a convenient 



approximation to the parameter family Gg that we will use in the proof of the main result Theorem 2.1 
From {H2) we have that there exists a positive constant Gi such that 

0<G,{x)<Gi, Va;e(0,l), Vee(0,eo)- (4.1) 

Moreover, for each x € [0, 1], we consider the function 

Go(a;) = minG(a;,?/) > 0. (4.2) 

We show that Go(-) is a continuous function in [0, 1]. Indeed, we will prove that 

|Go(a;)-Go(a;')| <sup|G(a;,y)-G(x',y)| Va;, x' G [0, 1]. (4.3) 



Consequently, the continuity of Go follows from the uniform continuity of G in y and inequality (4.3). 

Thus, let us prove (4.3). Given x and x' S [0,1], there exist y{x) and y{x') G K such that Go(a;) = 
G{x,y[x)) and Go(a;') — G{x' ,y{x'j). On the one hand, we have 

Go(a;) - Go(a;') = Go(a;) - G(a;, y{x')) + G(x, y{x')) - G(a;', y{x')) < G{x, y{x')) - G{x' , y{x')). (4.4) 

In a completely symmetric fashion we also obtain 

Go(a;') - Go(x) < G(x', y(a;)) - G(a;, y{x)). (4.5) 



Consequently, we obtain (4.3) from (4.4) and (4.5). 

Now, let us denote by the largest integer such that N^Le" < 1, where L is given in hypothesis (H2). 
Observe that N, ~ L-'^e'". Let 



Gne— rnin G ( a;, — ) 

xe[(Ti-l)Le°,nLe°] V e" / 



n = l,2...,iV, 



(4.6) 



and 7„_g € [(n — l)Le",nLe"] a point where the minimum (4.6 1 is attained, that is, G{'~fn,e, '^^) — G„^e 
where 7„ ^ does not need to be uniquely defined. By extension, let us denote by 70, e = and jN^+i.e ~ 1- 
Note that the set 

{70,e,7l,e, •••,7A'.+l,J (4-7) 
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defines a partition for tlie unit interval [0, 1]. Also, we have by definition that the segments 

{(7„,,,a;2) I G„,, < X2 < Gi}nn' = 0, 



for all n^l,2,...,N^. 

Consider also the step function 

f Gi,„ X e [0,71,,] 

Go(x) = <! IIia-}i{Gn,e,Gn+l,e}i 2: G [7„,e , 7„+i,e] , n = 1 , 2 . . . , iV, - 1 . (4.8) 

Lemma 4.1. We have 

l|Go - (5^11^00(0,1) -> ase->0. 



Proof. It follows from (i?2) and (4.3) that, for each 77 > 0, there exists eo > such that 

max{|G(x, y) - G(x', y)\, \Go{x) - Go(x')|} < V (4.9) 

whenever |a; — x'| < 2 Eq L and y e M. Now, for all x € [jn.e, 7n+i,e] we have 

Go{x) - Go(2:) = max{G„,,,G„+i,J - Go{x). 

Without loss of generality, we may assume Gl(x) = G„,e, that is, Gn.e > Gn+i.,,. Thus, 

Q(x)-Go(x) = G„,e-Go(a;) 

= G(7„,„7„,e/e")-Go(x) 

= G(7„,„7„,e/e")-Go(7n,.) + Go(7n,.)-Go(a;). (4.10) 
It follows from definition of Gq in ( |4.2[ ), that 

G(7n,o7n,e/e")-Go(7«,e)>0. 

Also, since G(a;, •) is Z(x)-periodic with |/(x)| < L, we have that there exist y{jn.e) G [0, ?(7n,e)] and k{jn.n) G 
N with yi-fn,e) + k{jn,e)lhn,e) € [(n - 1 ) i e" , n i e"] , such that 

Go(7n,e) = G(7„,e,y(7„,£)) = G{jn,e,yi,^n,e) + k{jn,e) li,^n,e))- (4.11) 

Consequently, we get from (4.6) and ( |4.11 1 that 

G(7„,., 7«,.A") - Go(7n,e) = G(7„,e,7n,e/e") - G((y + kl) e", (y + fcO e"/e") 

+G((2/ + k I) e", (y + A: e"/e") - G(7„., (y + fc e"/e") 
< G((y + A:0e",(y + fc0e7e")-G(7n,e,(2/ + fc0e"A") (4-12) 

since 

G(7n,e, 7n,e/e") " G((y + fc 6«, (y + e"/e«) < 0. 



Therefore, due to (4.10), (4.12) and (4.9), we obtain 

\Gt>{x) - Go{x)\ < G{{y + kl) e'',{y + kl) e'^/e") - G(7„,„ {y + kl) e"/e") + Go{ln,e) - Go{x) < 27j 

whenever x e [7„,e, 7„+i,e]- 

Then, since x G [0, 1] is arbitrary and U^li['yn,e,ln+i,e] = [0, 1], we conclude the proof. □ 

The following result will also be needed. 
Lemma 4.2. We have the following 

G4-)'^°77t/ G{;s)ds w*-L°"(0,l). (4.13) 



n-) Jo 
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Proof. We have to prove 

/ IgAx)- ['^^ G{x,s)ds\ip(x)dx^Q ase^O (4.14) 
Jo { Kx) Jo J 

for all ipe ii(0, 1). 

Since the set of step function is dense in i"'^(0, 1), and any step function is a linear combination of 
characteristic functions, we only need to show ( |4.14 1 for characteristic functions. Then, for < e < / < 1 
we consider the following characteristic function 



So, we have to estimate the integral 

/e,/= [{G^ix)- 



1 xe{ej) 
xiiej) ■ 



/ G(x, s) ds > dx 

Kx) Jo J 



as e > goes to zero. 

For this, let 77 > be a small number and let {e = a;o, a^i, a;„ = /} be a partition for the interval (e, /), 
and Xi be a fixed point in the interval Jj = [xi-i,Xi\, i — 1, such that 



sup sup \G{x,y) - G{xt,y)\ < -q. 



Observe that we can write 



where 



n „ 

Z2 J I ~ I(£Ty y G{xi,s)ds\dx 



i(xi) i.i(x) -\ 

G{x, s) ds > dx 



^''^^%lj\Wi)lo 



l-lix) J i-lix) 



G(x,s)ds——— / G(x,s)ds>dx. 

Kx) Jo J 

It is easy to estimate the integrals /f^, I* j: and ^ to obtain 

\llf\<v{f-e) 
\llf\<v{f~e) 

|/e'/l <Gi||^-;|U»(oa) (/-e) 

\l!j\<G, (l/L'^) \r-l\\L^^os){f~e) 



(4.15) 



where Gi, L and L' are the positive constants given by hypothesis {H2), and the function is the step 
function defined for each 77 > by 

l^{x) — l{xi) as Xi S Ji. 
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Since the inequalities (4.15) do not depend on e > 0, and — IW^o 



(0,1) 



as 77 — >■ uniformly in e, we 



have that j, j, ^ and ^ go to zero as 77 — > uniformly in e > 0. 

Hence, to conclude the proof, we just evaluate the integral y. But this is a straightforward application 
of the Average Theorem since Xi is a fixed point in Ji, and G{xi, ■) is a Z(a;i)-periodic function. Indeed, 



7^2 



E 



-1 J J, 



1=1 



G(i„x/e") 



K^i) Jo 



G(xi, s) ds > dx ^ as e ^ 0. 



□ 



4.2. Proof of Theorem 2.1, Here, we give a proof of the main result. Theorem 2.1 
Proof. The variational formulation of (2.2) is: find u*^ G iJ^(r2'^) such that 

I- 0x1 0x1 o'a;2 02:2 J Jqi 



(4.16) 



Taking = u*^ in expression (4.16) and using that ||/^||L2(f2e) < C, we easily obtain the a priori bounds 



dxi 



and - 



(9X2 



L2(n») 



< C. 



(4.17) 



In particular, we have 



<eC^Oase^O. 

0x2 L^(n') 

Let us observe that domain fJ*^ consists of two main parts. One of them is a highly oscillating domain 
and the other one is a non oscillating domain 571. To define these domains, we use the partition for the unit 
interval (0, 1) given by (4.7) and the step function Gg defined in (4.8). So, we consider the following open 
sets 

r!l = {(a:i,X2) G I XI e (0, 1), -6(xi) < X2 < Gg(xi)} 
= {(xi, X2) e I ^ (0, 1), Gg(xi) < X2 < a(xi)}. 

Notice that 

= Int (npjTF) . 



(4.18) 



Observe also that with the function Go defined by (4.2) we may also consider the set 
^0 = {(xi,X2) e I XI e (0, 1), -b{xi) <X2< Go(xi)}. 

which satisfies fin C fi'^. 



We want to pass to the limit in the variational formulation (4.16) for certain appropriately chosen test 
functions. In order to accomplish this, we rewrite it as follows 



f du"^ dip 1 dip \ 

1 ^ — ^ — — )dxidx2 

^ 0x1 0x1 ax2 0x2 J 



r du'^ dp 1 du"^ dp ) 

1 ^ — ^ — ^ — )dxidx2 

L 0x1 0x1 e'^ 0x2 0x2 J 



■ / u'^pdxidx2— / f'^ipdxidx2,'^p & H^{i^'^)- 



(4.19) 



Now, we pass to the limit in the different functions that form the integrands of (4.19). 

(a). Limit of u'^ in L^. 

It follows from (4.17) that u'^\ng € -ff^(f2o) and satisfies for all e > 



l"1lL2(f2o): 









and - 


dxi 


L2(no) e 


dx2 



< c. 
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Then, we can extract a subsequence of {w^lnol C H^{ilo), denoted again by u'^, such that 

u'^ ^ Uq w — H^{no) 
^ uo s- i?'*(rio) for all s e [0, 1) and 



dx2 



(4.20) 



s-L^{no) 



as e — > for some uq G iJ^(rio)- 



A consequence of the limits (4.20) is that Mo(a;i,a;2) does not depend on the variable X2- More precisely, 

duo , 



dx- 



-{xitX2) — a.e. 



(4.21) 



Also, due to (4.20), we have that the restriction of to the coordinate axis xi converges to mq. That is, 

for all s e [0, 1/2) where T = {(a;i, 0) € | xi e (0, 1)}. Consequently, we obtain 

— wo||l2(p) — > as e -> 0. 



Now, we can see that (4.22) implies the L^-convergence of to uq, that is 

- uo||l2(q.) -> as e -> 0. 



(4.22) 
(4.23) 



In fact, it follows from (4.22) that 



\\u''{xi,0) - UQ{xi)\\l2m,) = lu^i^xi^O) - uo{xi)f dx2dxi 

Jo J-Hxi) 

< C{G,b)\\u'~uo\\LHr) 
-> as e ^ 

where C(G, b) is a generic constant that depends on the functions G and b. Also, 



u''{xi,X2) - u''{xi,0) ^ - — {xi,s)ds 

Jo OX2 



and with Holder inequality. 



\u\xr,X2)-u'[x^M^ < 



du'- 
dx2 



{xi,s) 



ds \X2\ 



1 r'^eixi) 



Hence, integrating in fl*^ and using (|4.17|) to get 
||u'(a;i,a;2) -u'(a;i, 0)11^2(0.) 



1 i'Gc(xi) / rX2 



< 

lO J-b{xi) \Jq 

Therefore 



J-b{xi) 
2 \ 



\u'^{xi, X2) — u'^{xi, 0)1^ dxidx2 



dx- 



■(xi,s) 



ds \ \x2\ dx2dxi < C{G,b) 



dx2 



< eC{G,b) ^ as 0. 



L2(n«) 



- ^011^2(^2.) < \\u''{xi,X2) - u''{xi,0)\\L2(n') + \\u^ixi,0) - uo(a;i)|li2(s;2c) 
as e 0. 



(b). Limit of /^ 

Since ||/'|jL2(-j7e) < C independent of e, we have the function defined by 

r{xi) = / r{xi,x2)dx2 

J~b{xi) 



(4.24) 
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belongs to L^{0,1) and satisfies ||/'^||l2(o,i) ^ C{G,h) for some constant C{G,b) independent of e. Hence, 
via subsequences, we have the existence of a function / = f{xi) G L^(0, 1) such that 

r^f w-L\0,l). (4.25) 

Remark 4.3. Observe that in the case where f'^{xi,X2) = f{xi) then 

rixi) = (G(xi, xi/e") + b{xi)) /(xi) - p{xi) f{xi) w* - L°°(0, 1) 

where the function p is given by 



Pi^) = TT^ G{x,s)ds + bix), 

Jo 



which is the weak *-L°°(0, 1) limit of Ge{x) obtained in Consequently, we have that 

f{x)^p{x)f{x) xe(o,i). 



(4.26) 



(c). Test functions. 

Here, we define suitable test functions that will allow us to pass the limit in the variational formulation 



(4.19). For this, we use the definition of the open sets rii and £71 given in (4.18). 



n= 1,2,. 



For each (f) € H^{Q, 1) and e > 0, we define the following test functions in i7^(i7^) 
where is the rectangle (see Figure 4) 



x^Xxi,x2), {xi,x2) en%nQ^„ 

0(a;i), ixi,X2)eil.t 



Qn = {{X1,X2) |7n,£ < Xi < 7„+l,e, Gq{xi) < X2 < Gi} 

and the function is the solution of the problem 



dx\ 
dX" 



0, in 



0, on dQl\Tl 



, X\xi,X2) ^ (t){xi), on 
where F^ is the base of the rectangle, that is, 

= {{xi,Gl{xi)) : 7„,, <Xi< 7n+l,e}- 



(4.27) 



(4.28) 




Figure 4. Rectangle Q%^. 
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It follows from estimate (3.5 1 that 





1 


1 




dx\ 






dxl 



(4.29) 



Moreover, if we call = U^=iQn, we get = n rj% and we can define the function X'^ in fl'^ by 

X''{xi,X2) = X^{xi,X2) as {xi,X2) e n Q". 
Hence, we have by (4.29) that X'^ £ and satisfies the following inequality 



d^X' 



dx{ 



1 



d^X' 



dxl 



d^X' 



dx{ 



1 



Furthermore, we can show that 



Iv' - <t>\\L^(n') ^ as e ^ 0. 



(4.30) 



(4.31) 



We can argue as in (4.23). Indeed, since 



ip''{xi,X2) - (l>{xi) = ip''{xi,X2) - (/3'(xi,0) 



we have by \A.27\ and ( |4.30[ ) that 

,, ^ ^ , 2 



11^' - <C(G, 5) 



<C(G,6) 



9X2 



3a;2 
< C{G, b) e 



(xi, s) ds, 



L^iQ'nSl') 



as e ^ 0. 



(d). Passing to the limit in the weak formulation. 

Now, we pass to the limit in the variational formulation of the problem using the test functions ip"^ defined 



above. For this, we analyze the different functions that form the integrands in (4.19). 
• First integrand: 



\ ^ — ^ 1" — — \dxxdx2 as e ^- 0. 

L ax\ 0x1 0x2 0x2 J 

Indeed, it follows from ( |4.30 1 and a > that 

f du" dip" I du'' dip'^ , , f ( du" dX" 1 du" dX" \ , , 

\^ — ^ '^^^ — — \dxidx2 = / 7, 1-^^ — -T. — \dx1dx2 

L 0x1 0x1 0x2 0x2 J Jq,^ 0x1 0x1 0x2 0x2 J 



(4.32) 



1/2 



0. 



Second integrand: 

r du'^ d(p'^ 1 du"^ dip'^ 
I dxi dxi dx2 8x2 



^dxidx2 



To prove this, observe that using (4.27), we obtain 

dip' 



dxi 



dxi 



{Go{xi) + b{xi)) Uq{xi) (j)'{xi) dxi as e — >■ 0. 



1/2 



(4.33) 



(4.34) 



and 



dx2 



dx2 







for all e > 0. Hence, we have that 



[ ( du" dip" 1 du^dip"} , , f du" , x ,// x , 

/ i 7^ — -7. 1-^^ — -7 — \dxidX2= / -7 — {XI,X2)(P {xi)dxidx2 

Jq, L dxi dxi dx2 0x2 J Jw 0x1 
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dxi 



{xi,X2) 0'(xi) dxidx2 



dxi 



{xi,X2) (tJ^xi) dxidx2 



ni\ao 



dxi 



{xi,X2) (ty'ixi) dxidx2 



(4.35) 



Due to (4.20), we can pass to the limit as e — ^ in the first integral of the right side of (4.35) to 
obtain 



/no 

Also, we have that 



{xi,X2) dxidx2 



u'q{xi) (j)'{xi) dxidx2. 



no 



u'q{xi) (f)' (xi) dxidx2 



"o(^i) 'l>'{^i) 



Goixi) 



dx2 dxi 



b{xi) 



{Go{xi) + b{xi)) Uq{xi) (f>'ixi) dxi. 



(4.36) 



Now, we will get (4.34) if we prove that the remaining integrals of (4.35) tend to zero as e — > 0. 



We evaluate one of them, the other calculus is similar. 
From ( |4.17[ ), ( |4.18| and Remark [XT] we have that 



ni\f2o 



dxi 



{xi,X2) (f)' (xi) dxidx2 < 



dxi 



< C \ cb'ixif Goixi) - Glixi) 



dxi 



1/2 



< G \\(t>'\\L^{o,i) \\Go ~ G'o||^''i(Q 
as e -H- 0. 



(4.37) 



Therefore, we obtain (4.34) from (4.36) and (4.37). 



Third integrand: \i p(x) is defined in (4.26) then. 



Lp^ dxidx2 ~> I p(xi) uo(xi) 0(xi) dxi as e — > 



(4.38) 



To prove (4.38), observe that 

(fi"^ dxidx2 = I {u^ — uq) ip'^ dxidx2 + I uq {(p'^ — (/)) dxidx2+ I uq (p dxidx2 



From (4.23) and (4.31), we have 



/ (m*^ — Up) (p'^ dxidx2 — > and / ug {'p'^ — 4>) dxidx2 — > 0, as e — > 



In 

Hence, since 



uoixi) (j){xi) dxidx2 ^ / uo{xi) (l){xi) {Ge{xi) + b{xi)) dxi 



we get ( |4.38| ) from ( |4.13 ). 
Fourth integrand: 



(p'^ dxidx2 / /(xi) (/)(xi) dxi as e — 0. 



(4.39) 
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For this, let / e i^(0, 1) be the function defined in (4.24 1. Since 

/ if'^ dxidx2 — / {ip'^ — (j)) dxidx2 + / 4>dxidx2 
Jn' Jo.' Jsi' 



and 



(f) dxidx2 = 



b{xi) 



f''{xi,X2)dx2] (p{xi)dxi ^ / f{xi)(j){xi)dx 



we get ( |4.39| ) from (|2^, ( |4.25[ ) and ( |4.31[ ). 
Therefore, ifrom (4.32), ( |4.34 1, (4.38) and ( |4.39 1 we obtain the following limit variational formulation 

/•I . 

{(Go(xi) + 6(2:1)) +p(xi)uo(a;i)(/)(a;i)}dxi = / f{xi)(j){xi)dxi V0ei?^(O,l). 



(4.40) 

Since this problem is well posed, we obtain that the sequence {u'^}e>o is convergent and converges to the 
unique solution uq of (4.40). Thus we conclude the proof of Theorem 2.1 □ 

5. COMB-LIKE THIN DOMAINS 



We consider now. Type II thin domains as described in Section [2] and provide a proof of Theorem 2.5 



Proof. We will proceed as in the previous section to show this result. We will choose appropriate test 
functions to pass to the limit in the variational formulation of problem (2.8 1 that we rewrite it here as: find 
u" G ffi(rj') such that 

du'^ dip 
8x2 dx2 ■ 



f (du'^ dip 1 du'^ ^'^ \d d f / '^'^ ^ d 
Jq^ I dxi dxi dx2 dx2 / ^ ^ I dxi dxi dx-} dx-? i ^ ^ 

+ / u'^ipdxidx2— / f'^ipdxidx2,yip ^ H^{fl'^). 



(5.1) 



Again, as in the previous case, taking p = u'^ in expression (5.11 and using that ||/'^||L2(r2'=) ^ C"; '^s easily 
obtain the apriori bounds 



In particular, we have 



dW 









and - 


dxi 


L^{n') e 


3X2 



8X2 



L2(n=) 



L2(n«) 
<eC-J>Oase-^0. 



< C. 



(5.2) 



We extract a subsequence of {u'^|n_} C i?^(il_), denoted again by m*^, such that 

u" ^ uo s~ i?''(rj_) for ah s e [0, 1) and 



8X2 



(5.3) 



^0 s-L^{n^) 



as e — ^ for some uq E H^{il^). 



As in (4.211, it follows from (5.31 that uo{xi,X2) does not depend on the variable X2 and belongs to 

8uo 



H^{0, 1). Indeed, we can show that 



8x2 



(xi, 2:2) = a.e. 



(a). Limit of in L'^ifl''). 

First, we obtain the L^-convergence of u'^ to uq. More precisely, we show 

- uo||l2(o«) as e -> 0. 



(5.4) 
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For this, we assume without loss of generahty that 

C {{XI,X2) e M2 I Xi e (0, 1), 0<X2< b{xi)} 

and we define by 'symmetrization' the foUowing function u'^ in 51*^ by 



u''{Xi,X2) 



U''{xi,-X2), {xi,X2)eQ% 
U''{X1,X2), {xi,X2)efl-. 



(5.5) 



Consequently, it follows from (5.3 1 that 



Wu" - uo\\L2(n>) ^- as e ^ 0, 



and from (5.2), we have 



Let us denote by 



dxi 



and - 



3X2 



< C. 



(5.6) 



vf^ in ri'^. It is easy to see that w"^ = in 51- and w'^ satisfies 

lk'll_ffi(n;) = Wu" - -""11^1(12;) < Ci 

+ ) < C2 

Now let us show that ||w'^||i2(t2«) — as e — > 0, that is, ||w'^||L2(n^) as e -> 0. Suppose this is not 





2 


1 




dxi 






8X2 



(5.7) 



true and assume that \\w'^ 



This implies that the first eigenvalue of the problem 



c 1 > Co > at least for a subsequence e — > 0. Then we have that 





1 




2 


dxi 




dx2 





< 



C2 

Co 



C. 



1 92^^ 



dv"^ J 1 Sw'^ g 
dxi ^ 9a;2 



in rj^ 



on dn^XT 
w'(xi,0) = onT 



(5.8) 



satisfies Ae(i7i) < C, since J is the associated Raleigh quotient and F C Sfi^ is a nonempty open subset. 



But observe that fi^ = U^lj^fi^ ^ where all 17^ ^ are disjoint and identical, except for translations. Then, 
we can conclude X^(^l^) = \e{^n + ) ^■ 

Performing in _|_ the change of variables that transforms it into the fixed domain Qqj that is, (xi, X2) — > 
(xi/e" — nL, X2), we will have that Ae(ri^ _|_) is the first eigenvalue of the problem 



1 d'^v' 



1 d^v'^ 
1 dv^ 



+ v'^ ~ v'^ in Qo 



and therefore, 

Ae(J7^) = min ■ 



Qo 



9x1 2 



w'(xi,0)=0 onFo 



(5.9) 



dxidx2 



L\v'\ dxidx2 



eH\Qo),v'\r„^0} <C. 
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But this is impossible since Ae(17i) > -^ei{Qo) for a > 1 and e G (0, 1) where 



ei(Qo) = min • 







2 






( 




+ 

















dxidx2 



Iqo I"' I dxidx2 



is the first eigenvalue of the Laplace operator in Qq with homogeneous Dirichlet boundary condition in Fq 
and Neumann everywhere else. This eigenavalue is strictly positive by hypothesis (HQ). Thus we obtain 
@. 

(b). Test functions. 



(5.10) 



The test functions we are going to construct to pass the limit in the variational formulation (5.1 ) are very 
similar to the ones we constructed in Type I thin domains. Take (p G ^^^(0, 1), e > and define the following 
functions in 7J^(f2'^): 

( X^„{x,,X2), {x,,X2)enXnQ'„, n = l,2,...,N^^j^ 

^ 1 0(xi), {x,,x2)en^ 

where is the rectangle 

Ql = (nLe",(n + l)Le") x (0, G) 
and the function is the solution of the problem 

d^x^ 1 d^x^ . 

0, m 



dxi dx2 
BX^ 

— - 0, on dQ^^\K 
, X^{xi,X2) = (p{xi), on 
where is the base of the rectangle, that is, 

= (nLe", (n + l)Le") n 9Qo- 
As we showed in the previous section, we have using Lemma |3.H 

,a-l II ^'11 2 



(5.11) 



9V 


1 


1 


(9 V' 


dx\ 






dxj 



< Ce° 



L2(0,l) 



which implies that 

ll^1lL2(a^) + 
Moreover, we can show that 



(9 V 


2 


1 


9V 


dx\ 






dxl 



< c. 



- ^LHQ') ^ as e ^ 0. 



(5.12) 

(5.13) 
(5.14) 



We can argue as in (5.4 1. If it were not true, then there will exists a Cq > and a sequence (that we still 
denote it by e) such that \\ip'^ — </>|1l2(o«) > cq. But then, if we define w'^ — ip'^ ~ <f), we will have that 





2 


1 




2 


dxi 






8X2 





Co 



but with the same steps as we did in (a) this will contradict the fact that ei{Qo) > 0. 
(c). Pass to the limit. 
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Now we can pass to the limit in the variational formulation (5.1 1. First, we note that the convergences of 

(5.15) 



1 ^ — ^ ^ — ^ — \dxidx2 

L oxi oxi 0x2 0x2 J 



for e ^ follows from (5.121 and can be obtained as in (4.32). 
Also, from (5.3) and since Lp'^ = (j) \n we easily get 

r du'^ d^f 1 du'^ dtp'^ ~| 

\ ^ — ^ ^ — a — }dxidX2 

I 0x1 0x1 0x2 0x2 J 

Let us consider now the following technical result. 

Lemma 5.1. We have 

_ \Qo\ 



b{xi) Uq{xi) (/)' (xi) dxi as 



(5.16) 



/ dx2 q{xi) = 



L 



+ b{xi) w*-L°°{0,l). 



Proof. If we denote by x the characteristic function of the measurable open set Qq, extended periodically 
with respect to the first variable, we have by the Average Theorem that 



dX2 



xixi/e°',X2) dx2 



dxo 



IQol 

L 



\1j j ^(^'^2) dsj + &(xi) 
+ b{xi) Va;ie(0,l). 



u;* -i°°(0,l) 



□ 



Moreover, 



u*^ (/j*^ dxidx2 



(u*^ — uq) ip'^ dxidx2 + I Uq {if'^ — (j)) dxidx2 + / UQ(j)dxidx2- 

^llL2(f2«) ~^ 0. The last integral 



and the first two integrasl go to since — u'^\\L^{n^) — > and \\(p' 
satisfies, 



uq{xi) 4>{xi) dxidx2 



Uq{xi) (l){xi) 



dx2 dxi, 



S'{xi) 



q{xi)uo{xi) (p{xi)dxi 



where we have used Lemma l5.1 
Finally, we have 



Lp'^ dx\dx2 



dxidx2 + / f'^4>dxidx2 



but the first integral goes to 0. Moreover, with the hypothesis of the theorem, we get for the second integral 



(j)dxidx2 = 



0= 







Therefore, we obtain from the estimates above that 
{b{xi) u'q[xi) (f>'{xi) + q{xi) uo{xi) 0(a;i)} dxi 





r{xi,X2)dx2] (l){xi)dxi / f{xi)(j){xi)dxi 



f{xi)<l>{xi)dxi V0ei/i(O,l). 



(5.17) 



Since this problem has a unique solution, then we obtain that the sequence {u'^}e>o is convergent and 
converges to the unique solution uq of (5.171. 

□ 
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